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Abstract  We  consider  the  media-access  control  problem  for 
nodes  with  heavy  traffic  in  single-destination  wireless  net¬ 
works.  We  assume  that  each  source  transmits  in  each  time 
slot  according  to  a  transmission  probability,  which  is  a  con¬ 
tinuous  value  between  0  and  1.  Our  goal  is  to  determine 
the  values  of  the  transmission  probabilities  so  that  the  net¬ 
work  throughput  is  maximized.  In  this  paper,  we  show  that 
the  maximum  throughput  is  achieved  only  if  these  values 
are  either  0  or  1.  We  obtain  closed-form  results  for  optimal 
throughput  for  networks  that  operate  under  a  homogenous 
situation  in  which  the  expected  value  of  the  received  power  at 
the  destination  is  the  same  for  each  source.  We  then  extend 
our  studies  to  more  general  networks,  which  rely  on  exhaus¬ 
tive  search  for  the  optimal  set  of  transmissions.  The  search 
has  exponential  complexity  and  is  feasible  only  for  networks 
with  small  or  moderate  sizes.  Thus,  we  also  develop  heuristic 
algorithms,  which  have  polynomial-time  complexity  and  are 
suitable  for  large  and  general  networks. 

1  Introduction 

We  study  the  media-access  control  (MAC)  problem  for  one 
of  the  most  basic  forms  of  wireless  networking:  a  network 
with  K  sources  and  a  single  destination.  An  example  net¬ 
work  is  shown  in  Fig.  1.  Such  a  single-destination  may  be 
part  of  a  much  larger  multiple-source  multiple-destination 
network  (e.g.,  a  cellular  network  or  a  sensor  network).  The 
network  operates  in  a  complex  environment  that  includes 
receiver  noise,  other-user  interference,  and  channel  fading. 
Assume  that  the  receiver  has  multi-packet  reception  (MPR) 
capability,  i.e. ,  one  or  more  packets  may  be  decoded  suc¬ 
cessfully  if  the  signal-to-interference-plus-noise  ratio  (SINR) 
exceeds  a  specified  threshold  [3,  4], 

We  focus  on  the  MAC  problem  for  the  uplink  channel 
(i.e.,  the  channel  from  the  sources  to  the  destination),  and 
we  are  interested  in  how  each  of  the  K  sources  transmits 
in  each  time  slot.  Here,  we  study  a  model  that  is  different 
from  existing  models  such  as  random  access  (e.g.,  Aloha), 
CSMA/CD,  CSMA/CA  (e.g.,  802.11)  and  CDMA  [2,  6,  7, 
8].  In  particular,  our  model  assumes  that  source  i  transmits 
in  each  time  slot  according  to  a  transmission  probability  pi , 
which  is  a  continuous  value  between  0  and  1.  Our  goal  is 
to  determine  the  vector  (pi,P2,  ■  ■  ■  ,Pk)  so  that  the  network 
throughput  is  maximized. 

In  this  paper,  we  show  that  the  maximum  throughput  is 
achieved  only  if  the  transmission  probabilities  are  either  0 
or  1.  This  result,  which  holds  for  any  form  of  fading  and 
receiver  noise,  implies  that  the  search  space  for  the  optimal 
probabilities  is  reduced  from  an  infinite  and  continuous  space 
to  a  finite  set  of  points.  Thus,  our  problem  is  equivalent 
to  that  of  determining  a  subset  of  transmission  nodes  that 
maximizes  the  throughput. 


We  obtain  the  closed-form  optimal  solution  for  networks 
that  operate  under  a  homogenous  situation  in  which  the  ex¬ 
pected  value  of  the  received  power  at  the  destination  is  the 
same  for  each  source.  This  case  may  arise  when  the  sources 
are  located  at  the  same  distance  from  the  destination,  or  al¬ 
ternatively  when  power  control  is  used.  We  then  extend  our 
studies  to  more  general  networks,  which  rely  on  exhaustive 
search  for  the  optimal  set  of  transmission  nodes.  The  search 
has  exponential  complexity  and  is  feasible  only  for  networks 
with  small  or  moderate  sizes.  Thus,  we  also  develop  heuristic 
algorithms,  which  have  polynomial-time  complexity  and  are 
suitable  for  large  and  general  networks.  Additionally,  our 
numerical  results  show  that  these  heuristics  perform  almost 
identically  to  the  optimal  solutions  obtained  by  expensive 
exhaustive  search. 


Fig.  1  A  network  with  K  sources  (Si)  and  a  single  destination  ( D)\ 
Source  «S)  transmits  with  probability  pi  in  a  time  slot. 

2  Network  Model  and  Assumptions 

We  consider  a  stationary  wireless  network  that  has  K 
sources,  denoted  by  Si,  S2, ■■■ ,  Sk,  that  transmit  their  traf¬ 
fic  to  a  common  destination,  denoted  by  D.  An  example 
network  with  K  =  6  sources  is  shown  in  Fig.  1.  We  assume 
the  following: 

•  The  nodes,  whose  locations  are  known  and  fixed,  are 
equipped  with  omnidirectional  antennas. 

•  The  destination  can  successfully  receive  more  than  one 
transmission  at  a  time,  i.e.,  it  has  MPR  capability. 

•  Each  source  can  communicate  directly  with  the  destina¬ 
tion.  Routing  is  not  considered  in  this  paper. 

•  The  traffic  is  heavy  in  the  sense  that  each  source  always 
has  traffic  to  transmit,  i.e.,  its  transmission  queue  is  never 
empty. 

•  Time  is  divided  into  slots.  The  traffic  is  expressed  in 
terms  of  fixed-size  packets  such  that  it  takes  one  time 
slot  to  transmit  one  packet. 

•  Our  primary  performance  measure  is  sum  throughput, 
which  is  the  average  number  of  packets  that  are  success¬ 
fully  received  by  the  destination  per  time  slot.  We  do  not 
address  issues  such  as  time  delays  and  stability  analysis 
in  this  paper. 

•  Source  Si  transmits  with  probability  pi  in  each  time 
slot,  0  <  Pi  <  1.  Our  goal  is  to  determine  the  vector 
(pi,P2, . . .  ,pk)  that  maximizes  the  throughput. 
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We  assume  that  a  packet  is  successfully  received,  even  in 
the  presence  of  interference  and  noise,  as  long  as  its  SINR 
exceeds  a  given  threshold  [3,  4].  More  precisely,  suppose  that 
we  are  given  a  set  H  of  sources  that  transmit  in  the  same 
time  slot,  and  S  G  H.  Let  Prx(S,  D)  be  the  signal  power 
received  by  node  D  for  node  S,  and  let  SINR(5,  D)  be  the 
SINR  at  D  for  the  transmission  from  S,  i.e. , 


SINR  (S',  D) 


_ fixCS,  £>) _ 

-Pnoise(-D)  +  PIX(U,D) 

UGH\{S} 


where  Pnoise(D)  denotes  the  receiver  noise  power  at  D.  We 
assume  that  a  packet  transmitted  by  S  is  successfully  re¬ 
ceived  by  D  if 

SINR(S,  D)  >  /3  (1) 

where  /3  >  0  is  a  threshold  at  D ,  which  is  determined  by 
application  requirements  and  the  properties  of  the  network. 
When  /?  <  1  (e.g.,  in  spread-spectrum  networks),  it  is  pos¬ 
sible  for  two  or  more  transmissions  to  satisfy  (1)  simultane¬ 
ously. 

The  wireless  channel  is  affected  by  fading,  as  described 
below.  Let  Ptx{S)  be  the  transmit  power  at  node  S,  and 
r(S,  D)  be  the  distance  between  nodes  S  and  D.  When  node 
S  transmits,  the  power  received  by  node  D  is  modeled  by 

PIX(S,D)  =  A(S,D)g(S,D) 

where  A(S,  D)  is  a  random  variable  that  incorporates  the 
channel  fading.  We  refer  to  g(S,  D)  as  the  “received  power 
factor,”  which  depends  on  r(S,D)  and  Ptx(S).  For  far-field 
communication  (i.e.,  when  r(S,D)  1),  we  have 


is  computed  as  follows.  To  simplify  the  notation,  we  define 
Qi  =  1  -  Pi- 

Let  H  be  a  set  of  sources,  and  let  Ch{Su  D)  be  the  prob¬ 
ability  that  a  packet  from  source  5)  is  successfully  received 
by  destination  D  in  a  time  slot,  given  that  all  the  sources  in 
H  simultaneously  transmit  in  the  slot. 

First,  consider  the  simplest  case  of  K  =  2  sources.  The 
throughput  is  T(pi,p2)  =  T1+T2,  where  T)  is  the  throughput 
contributed  by  ,S{ ,  i  =  1,2.  The  term  T)  is  the  sum  of  2 
components:  The  1st  component  is  the  probability  that  a 
packet  from  source  Si  is  successfully  received  by  destination 
I?  in  a  time  slot,  given  that  both  S 1  and  S 2  transmit  in  the 
time  slot  (which  occurs  with  probability  pip2).  Thus,  the  1st 
component  is  Cislts2}(Si,  D)p3p2.  The  2nd  component  is 
the  probability  that  a  packet  from  Si  is  successfully  received 
by  D  in  a  time  slot,  given  that  only  Si  transmits  in  the 
time  slot  (which  occurs  with  probability  Piq2).  Thus,  the 
2nd  component  is  (Si,  D)piq\.  We  then  have  Tf  = 
c{SuS2}(Si,  d)PiP2  +  C{Sl}{Si,D)piq2.  Similarly,  it  can  be 
shown  that  T2  =  C'{s1,s2}(-S'2,  D)PiP2+C{s2}  (S2,  D)qip2.  To 
summarize,  the  throughput  for  the  case  of  K  —  2  is 


T(pi,p2)  =  C{Sl,s2}(Si,D)pip2  +C{Sl}(Si,D)piq2 
+  C{Si,S2}(‘SI2,  D)pip2  +  C{S2}(S2,  D)qip2 

It  can  be  verified  that  the  right  hand  side  of  (4)  equals 
E(L)  =  Pr {L  =  1}  +  2Pr{L  =  2},  where  L  G  {0,1,2}  is 
the  random  variable  representing  the  number  of  successful 
transmissions  in  a  time  slot.  Thus,  the  throughput  T(pi,p2) 
computed  in  (4)  is  indeed  the  average  number  of  successful 
transmissions  in  a  time  slot. 

For  the  case  of  K  =  3  sources,  it  can  be  shown  that  the 
throughput  is  given  by 


g(S,D)  =  Ptx(S)r(S,D)~a  (2)  T  (pi,p2,Ps)  = 


where  a  is  the  path-loss  exponent  whose  typical  values  are 
between  2  and  4.  A  simple  approximate  model  for  both  near¬ 
field  (i.e.,  when  r(S,D)  <  1)  and  far-held  communication  is 

g(S:D)  =  Ptx(S)[r(S:D)  +  l]-a  (3) 

where  the  expression  r(S,  D)  +  1  is  used  to  ensure  that 
g(S,D)  <  Ptx(S).  Under  Rayleigh  fading,  A(S,D)  is  ex¬ 
ponentially  distributed  [9,  p.  36]. 

Our  goal  is  to  study  methods  for  accomplishing  the  com¬ 
munication  between  the  sources  and  the  destination,  and  to 
evaluate  the  resulting  performance.  Existing  methods  such 
as  TDMA,  802.11,  and  Aloha  do  not  allow  simultaneous 
transmissions,  i.e.,  there  is  exactly  one  transmission  and  no 
other- user  interference  in  each  time  slot.  In  this  paper  we 
consider  MPR  approaches,  as  described  in  the  following  sec¬ 
tions,  under  which  more  than  one  transmission  is  allowed  in 
a  time  slot. 

3  Throughput  Evaluation 

Recall  that  source  Si  transmits  with  probability  pi  in  each 
time  slot,  0  <  Pi  <  1.  Our  goal  is  to  determine  the  vec¬ 
tor  {pi,p2,  ■  ■  ■  ,Pk)  that  maximizes  the  network  throughput, 
which  is  the  average  number  of  packets  that  are  successfully 
received  by  the  destination  per  time  slot.  The  throughput 


C {s^iSi,  D)piq2qz  +  C{s2}(S2,  D)qip2q3  +  C{s3}(S3,  D)q1q2p3 
+  C{Si,s2}(Si>  D)pip2q3  +  C{Slis2}(^2,  D)pip2q3 
+  C {s1,s3}{S\1  D)piq2p3  +  C{s1,s3}(^3,  D)piq2p3 
+  C{S2,S3}(‘5,2,  D)qip2p3  +  C{s2ts3}(S3,  D)qip2p3 
+  C{Si, s2, s3}(Si,  D)pip2p3  +  C{Sl,s2,s3}{S 2,  D)pip2p3 

+  C{s1,s2,s3}[S3,D)pip2p3 

For  the  general  case  of  I\  sources,  it  can  be  shown  that 
the  throughput  is  given  by 

T(pi,P2,...,Pk)  =  n  Pi  n  yi  Ch(S,D) 

H^toSi^H  Sj£H  SGH 

(5) 

with  the  convention  Ils  eaU  —  Pj)  =  1- 

Direct  computation  of  the  throughput  T(pi,p2, . . .  ,pk) 
as  given  in  (5)  requires  the  summation  of  terms  whose  num¬ 
ber  increases  rapidly  as  K  increases.  As  seen  above,  the 
number  of  terms  grows  from  4  for  K  =  2  to  12  for  K  =  3.  It 
can  be  shown  in  general  that  the  number  of  terms  is  K 2K~1. 
In  the  following,  we  show  that  at  the  point  of  maximized 
throughput,  pi  must  be  0  or  1  for  all  i,  i.e.,  the  computation 
of  throughput  can  then  be  greatly  simplified. 
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4  Reduced  Search  Space  for  Optimal  Solution 

The  following  theorem  is  general,  because  it  holds  for  ar¬ 
bitrary  values  of  C#  (5),  .D)  used  in  the  evaluation  of  the 
throughput  T(pi,p2,  ■  . .  ,pk)  as  given  in  (5),  i.e. ,  the  theo¬ 
rem  does  not  depend  on  the  form  of  receiver  noise,  power 
attenuation,  and  fading. 

Theorem  1  Assume  that  each  source  Si  transmits  with 
probability  pi  in  a  time  slot,  0  <  Pi  <  1.  The  throughput 
T(pi,P2,  ■  ■  ■  ,Pk)  is  then  maximized  for  some  binary  vector 
(pi,P2,  •  •  •  ,Pi<),  i.e. ,  pi  must  be  0  or  1,  for  all  i.  Thus,  the 
optimal  solution  is  confined  to  a  finite  set  of  2K  elements. 

Proof  Consider  an  arbitrary  source  Si,  which  trans¬ 
mits  with  probability  pi  in  a  time  slot,  0  <  p,;  <  1 . 
For  each  (A'  -  l)-tuple  (pi,p2, . . .  ,p»-i,p»+i,  •  •  • ,Pk),  let 
us  study  the  dependence  of  the  throughput  7)  (pl )  = 
T(px,p2,  .  ■  ■ , Pi-i, Pi, Pi+i, . . .  ,pK)  on  pi.  First,  note  that 
Ti{pi)  is  a  weighted  sum  of  terms,  each  of  which  con¬ 
tains  either  p,  or  (1  —  pi).  Thus,  for  each  ( K  —  l)-tuple 
{pi,P2,  ■  ■  ■ , Pi-i, Pi+i, ,  ■ .  ,pk),  Ti(pi)  varies  according  to  a 
straight  line. 

Because  0  <  pi  <  1,  we  see  that  this  straight  line  starts 
at  pi  =  0  and  ends  at  pi  =  1.  Note  that  the  slope  of  this 
straight  line  depends  on  the  values  of  pj  and  Cn{Sk,D)  for 
all  j  i  and  all  sources  Sk  ■  Because  Ti(pi)  varies  according 
to  a  straight  line,  its  maximum  must  occur  at  one  of  its 
endpoints,  i.e.,  T)(pj)  <  max{X((0), Ti(l)}.  Thus,  we  have 


T (pi , P2 ,  ■  ■  ■ , Pi-i, Pi, Pi+i,  ■  ■  ■  ,Pk) 

<  max{T(pi,p2,  •  •  • , Pi-i,  0,Pi+i,  ■  ■  ■ ,Pk ), 

T(pi,p2, . . . , Pi-1,  l,Pi+ii  ■  ■  ■ ,PK )} 


(6) 


Reasoning  as  above,  the  terms  on  the  RHS  of  (6)  can  be 
further  expanded  as 
T{pi,P2,  ■  ■  ■ , Pi-i, 0,pi+i,  ■  ■  ■  ,pK) 


<  max{T(pi,p2,  •  ■  •  ,0,0,pi+i, . . .  ,pK), 

T(pi,P2,  ■  ■  ■  ,l,0,pi+i, . . .  ,pK)} 


and 

T{pi,P2 


,  Pi —  1 ,  ^-,Pi+l,  •  •  •  ,PK ) 

<  max{T(pi,p2,  •  ■  •  ,0,  l,P;+i, . . .  ,pK), 

T(pi,p2, . . .  ,l,l,pi+i, . . .  ,pK)} 

From  (6),  (7),  and  (8),  we  then  have 

T(pi,p2, . . . , pi-i, Pi, Pi+i, . . .  ,pK) 

<  max{T(pi,p2,.  ■ .  ,0,0,pi+i, . . .  ,pK), 

T(pi,p2,  ■ .  ■  ,l,0,pi+i, . . .  ,pK), 
T(pi,p2,  ■ . .  ,0,l,pi+i, . . .  ,pK), 
T(pi,P2,  ■  ■  ■  ,l,l,pi+i,  ■  ■  ■  ,Pk)} 


(7) 


(8) 


(9) 


By  expanding  the  terms  on  the  RHS  of  (9)  repeatedly,  it 
can  be  shown  that 


T{pi,P2,  ■  ■  ■  ,Pk)  <  max{T(0, 0, . . . ,  0), . . .  ,T(1, 1, . . . ,  1)} 


i.e.,  T(pi,p2, . . .  ,pk )  is  maximized  at  one  of  the  binary  vec¬ 
tors  (0,0,  ...,0),...,(1,1,...,1). 

For  the  case  of  K  =  2,  the  proof  is  especially  simple, 
because  by  applying  (8)  repeatedly,  we  have 

T{pi,P2)  <  max{T(0,p2),7’(l,p2)} 

<  max{max{T(0, 0),  T(0, 1)},  max{T(l,  0),  T(l,  1)}} 
=  max{T(0, 0),  T(  0, 1),  T(l,  0),  T(l,  1)}  □ 

Note  that  pi  =  0  means  source  Si  never  transmits,  and 
Pi  =  1  means  source  S',;  always  transmits.  Theorem  1  states 
that,  when  the  throughput  is  maximized,  pi  must  be  either  0 
or  1  for  all  sources  S,;,  i.e.,  only  a  subset  of  sources  are  allowed 
to  transmit  and  all  others  are  blocked  from  transmissions. 
Thus,  the  maximized  throughput  can  be  simplified  to 

T=^Cff(S,D)  (10) 

sgh 

where  H  is  a  set  of  sources  that  transmit  in  the  time  slot 
(i.e.,  their  transmission  probabilities  are  1).  As  seen  later, 
Theorem  1  can  be  adapted  to  allow  all  the  sources  to  transmit 
in  multiple  time  slots  (see  Remarks  1  and  2). 

Theorem  1  is  valid  for  any  form  of  receiver  noise  and  chan¬ 
nel  fading.  For  the  rest  of  this  paper,  we  focus  on  the  case  of 
Rayleigh  fading,  for  which  the  following  result  (whose  proof 
is  given  in  [1,  5])  provides  the  exact  formula  for  Ch(S,D), 
which  depends  on  the  receiver  noise,  channel  fading,  receiver 
threshold,  and  other-user  interference. 


Theorem  2  Suppose  that  the  fading  between  a  transmitting 
node  S  and  a  receiving  node  D  is  modeled  as  a  Rayleigh 
random  variable  Yg  with  parameter  v(S,D).  For  S  ^  U, 
assume  that  Yg  and  Yjj  are  independent.  Given  that  all 
the  nodes  in  H  simultaneously  transmit  in  a  time  slot,  the 
probability  that  a  packet  from  S  is  successfully  received  by 
D  is 


Ch(S,D ) 


/  pPnoise(D)  \ 

V  V(S,  D)g(S,  D) ) 


n 


1  +  /3 


U£H\{S} 


v{U,D)g{U,D) 
V(S,  D)g(S,  D) 


where  f3  and  Pnoise(D)  are  the  required  SINR  threshold  and 
the  receiver  noise  power  at  D ,  respectively,  and  g(S,  D)  is 
the  received  power  factor  defined  in  (2)  and  (3). 


5  Networking  in  Homogenous  Situation 

We  now  assume  that  the  network  operates  in  a  homoge¬ 
nous  situation,  meaning  that  packets  from  all  sources  have 
equal  chance  of  being  received  successfully,  i.e.,  Cn{Si,  D )  = 
Ch(Si,D)  for  all  subsets  of  sources  H  and  all  i.  From  (10), 
the  optimal  throughput  is  simply  T  =  \H\Ch(Si,  D)  for 
some  subset  of  sources  H. 

For  the  case  of  Rayleigh  fading,  it  is  sufficient  to  assume 
that  v(Si,  D)g(Si,  D)  =  v(S\,  D)g(S±,  D)  for  all  i,  where 
v(Si,  D)  is  the  Rayleigh  fading  parameter  and  g(S,  D )  is  the 
received  power  factor,  as  defined  in  (2)  and  (3),  which  de¬ 
pends  on  the  distance  r(S,  D )  and  the  transmitted  power 
Ptx(S')-  This  case  may  arise  when  all  the  sources  transmit 
with  the  same  power  and  are  located  the  same  distance  from 
the  destination  (as  shown  in  Fig.  2),  or  alternatively  when 
power  control  is  used.  Thus,  the  expected  value  of  the  re¬ 
ceived  power  at  the  destination  is  the  same  for  each  source. 
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Theorem  3  Assume  that  v(Si,  D)g(Si,  D)  =  v(Si,  D)g(Si,  D) 
for  1  <  *  <  K,  and  also  assume  that  to  sources  simultane¬ 
ously  transmit  in  a  time  slot,  1  <  to  <  K.  The  throughput 
is  then  given  by 


exp  - 


T  =  m- 


/^-^noise  A 

v(S1,D)g(SllD)) 

(1  +  /3)m_1 


which  is  maximized  when 


m 


1 


1 

-ln(l  +  /?) 
K 


+  0.5J 


if  /3  >  1 

if  ex'K  -  1  <  £  <  1 
if  0  <  /3  <  ex!K  —  1 


Proof  Let  H  be  a  set  of  sources  that  transmit  in  a  time 
slot,  and  let  m  =  \H\.  Theorem  2  implies  that 


CH(Si,D)  =  CH(S1,D) 


avn  I _ 4  Pnoise _  | 

LXPy  v(S1,D)g(S1,D)  J 

(1  +  /3)m_1 


The  total  throughput  from  the  m  transmissions  is  then 


i.e. ,  the  optimum  is  to  =  1  for  (3  >  1.  □ 

Remark  1 

(1)  Theorem  3  implies  that  the  optimal  number  of  transmis¬ 
sions,  m ,  is  independent  of  Pnoise  for  all  /3  >  0.  Further,  for 
/3  >  1,  the  optimum  is  m  =  1,  independent  of  the  network 
size  K. 

(2)  Let  us  compare  the  throughput  performance  of  our 
method  with  the  well-known  slotted  Aloha  for  the  case 
of  Pnoise  =  0.  The  throughput  of  the  slotted  Aloha  is 
PAioha  =  1/e  =  0.368  for  sufficiently  large  K  [7].  Theorem  3 
implies  that  the  throughput  T  of  our  method  is  between  1 
and  K.  Note  that  TAioha  is  independent  of  (3,  while  T  is 
highly  dependent  on  /?.  Further,  the  slotted  Aloha  is  dis¬ 
tributed,  while  our  method  is  centralized. 

(3)  Theorem  3  specifies  the  optimal  number  of  sources,  m, 

that  transmit  in  each  time  slot.  For  the  case  of  m  <  K, 
we  can  use  a  round-robin  scheme  to  allow  all  the  sources  to 
transmit  in  multiple  time  slots  (i.e.,  a  method  for  fairness). 
For  example,  suppose  that  the  total  number  of  sources  is 
K  =  5,  and  the  optimal  number  of  sources  that  transmit  in 
each  time  slot  is  m  =  3.  All  the  sources  can  transmit  in  mul¬ 
tiple  time  slots  as  follows:  S\,  S2,  S3  (in  slot  1),  S4,  S5,  S  1  (in 
slot  2),  S2,S3,S4  (in  slot  3),  S5,Si,  S2  (in  slot  4),  S3,S4,S5 
(in  slot  5),  etc.  □ 


T(m) 


a-v-Tt  l _ /^Pnoise _  \ 

eXPl  v(Su D)g(Su D)  I 

TTl - - - — 

(1  +  !3)m~l 


By  considering  m  as  a  continuous  variable,  taking  the 
derivative  of  T(m),  and  then  setting  the  results  to  zero,  it 
can  be  shown  that  m  =  1/  ln(l  +  (3). 

Note  that  m  in  our  original  problem  is  an  integer  between 
1  and  I\.  Thus,  m  =  1  if  l/ln(l  +  /3)  <  1,  and  m  =  K  if 
l/ln(l  +  j3)  >  K.  Further,  m  «  l/hi(l  +  /3)  when  1  < 
l/ln(l  +  /3)  <  K.  In  this  case,  our  numerical  studies  suggest 
that  to  =  [1/  hi(l  +  /?)  +  0.5J .  To  summarize,  the  optimal  m 
is  given  by 


1 

[1/  ln(l  +  /3)  +  0.5J 

I< 


if  1/  ln(l  +  /3)  <  1 
if  1  <  1/  ln(l  +  j3)  <  K 
if  1/  ln(l  +  (3)  >  K 


which  is  equivalent  to 


1 


1 


-ln(l  +  (3) 
I\ 


if  /3  >  e  —  1 

if  ex!K  -l<p<e-l 
if  0  <  [3  <  ex!K  —  1 


Now  assume  that  (3  >  1.  It  can  be  shown  by  induction 
that  m  <  2m_1  for  m  >  1.  Using  this  inequality,  we  then 
have 

rp(-l\  (  (3P-n- 

T(l)  =  exp  -- 


>  TO- 


exp 


v(S1,D)g(S1,D) 

( _ -Pnoise  A 

V  v(.S1,D)g(S1,D) ) 


>  m- 


exp 


(- 


2m~1 

1 3Pn , 

v(S1,D)g(S_ 


1  ,D)) 


(1 +  PY 


=  T(m) 


Fig.  2  The  sources  lie  on  a  circle  centered  at  the  destination. 

6  Networking  in  Non-homogenous  Situation 

We  now  assume  that  the  network  operates  in  a  non- 
homogenous  situation,  meaning  that  packets  from  the 
sources  may  have  unequal  chance  of  being  received  success¬ 
fully,  i.e.,  CniSi,  D)  ^  Ch(Si,D)  for  some  set  of  sources  H 
and  some  i.  As  seen  in  Section  4,  Theorem  1  implies  that, 
given  a  general  wireless  network  of  K  sources  and  a  destina¬ 
tion,  our  problem  reduces  to  finding  the  optimal  subset  H* 
of  sources  that,  when  they  transmit  simultaneously  in  a  time 
slot,  maximizes  the  network  throughput 

T=  Y.  Ch*(S,D) 

sgh * 

For  a  given  set  of  K  sources,  there  are  2^  —  1  non-empty 
subsets  of  sources.  By  comparing  the  throughput  values 
among  these  2K  —  1  subsets,  the  optimal  subset  H*  of  sources 
can  be  found.  Thus,  exhaustive  search  has  0(2^)  complex¬ 
ity. 

In  the  following,  we  present  2  heuristic  algorithms  of 
greatly  reduced  complexity.  Our  numerical  evaluation  for  nu¬ 
merous  network  instances  of  small  and  moderate  size  shows 
that  they  produce  optimal  or  nearly  optimal  solutions.  In 
both  of  these  algorithms,  candidate  transmission  sets  are 
produced  by  choosing  an  initial  transmitting  node,  and  then 
adding  additional  nodes  into  the  time  slot  until  the  resulting 
throughput  does  not  increase  further.  The  best  such  trans¬ 
mission  set  is  the  one  that  results  in  the  largest  value  of 
throughput.  The  algorithms  for  constructing  the  candidate 
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transmission  sets  differ  in  terms  of  the  order  according  to 

which  the  sources  are  added  into  the  time  slot. 

Algorithm  1:  Systematic  Node  Selection 

•  Label  the  K  sources  according  to  decreasing  values  of 
their  throughput  (under  the  assumption  that  each  node 
transmits  alone  in  a  time  slot),  i.e. ,  Si  has  the  highest 
throughput  and  Sk  the  lowest.  This  can  be  done  via 
sorting,  which  has  O(AMogA')  complexity. 

•  For  each  i,  1  <  i  <  K,  form  the  following  candidate 
transmission  set  of  consecutive  sources 

A,;  —  {$£,  <S)-|-1,  •  •  .  ,  } 

where  i  +  rii  <  K.  Note  that  there  are  no  gaps  in  A;,  i.e., 
|Aj|  =  m  +  1.  The  sources  Si,  Si+i, . . . ,  Si+ni  are  added 
one  by  one  into  the  time  slot,  where  Si+ni  is  the  last 
source  beyond  which  the  throughput  does  not  increase 
further.  Let  T(Aj)  be  the  throughput  obtained  when  the 
sources  in  Aj  transmit  in  the  time  slot. 

•  Finally,  we  choose  the  transmission  set  A*  that  yields 
the  maximum  throughput,  i.e.,  the  chosen  A*  satisfies 
T(A*)  =  max{T(A,;),  1  <  *  <  A'}.  It  can  be  shown  that 
the  overall  complexity  of  Algorithm  1  is  0(AT2). 


Note  that  the  candidate  transmission  sets  for  Algorithm  1 
contain  no  gaps,  when  the  sources  are  ordered  from  highest 
expected  received  power  to  lowest.  Restricting  the  class  of 
policies  that  must  be  examined  in  this  manner  reduces  the 
complexity  of  the  algorithm. 


Algorithm  2:  Greedy  Node  Selection 

As  with  Algorithm  1,  K  candidate  transmission  sets  are 
constructed,  each  using  one  of  the  K  sources  as  a  starting 
point.  Each  additional  source  is  chosen  with  the  greedy  ob¬ 
jective  of  maximizing  the  throughput  in  the  slot.  Unlike 
Algorithm  1,  the  numbering  of  the  sources  in  this  algorithm 
is  arbitrary,  i.e.,  unrelated  to  their  throughput  in  isolation. 
The  basic  idea  is,  at  each  step,  to  add  the  source  that  yields 
the  maximum  throughput  into  the  time  slot.  More  specifi¬ 
cally, 

•  For  each  i,  1  <  i  <  K ,  perform  the  following  procedure: 
—  Step  1:  Add  source  S)  as  the  initial  source  into  the 

time  slot. 

—  Step  k:  Add  into  the  slot  the  source  that  yields  the 
maximum  throughput  computed  up  to  this  step. 

—  The  algorithm  stops  when  the  throughput  does  not 
increase  further. 

•  Repeat  the  above  procedure  with  S)  as  initial  source, 
for  each  i,  1  <  i  <  K.  Let  T(Si)  be  the  resulting 
throughput  obtained  when  S,  is  the  initial  source.  We 
then  choose  the  initial  source  that  yields  the  maximum 
throughput,  i.e.,  the  chosen  initial  source  S)*  satisfies 
T(Si »)  =  max{T(S'i),  1  <  i  <  K}.  It  can  be  shown  that 
the  overall  complexity  of  Algorithm  2  is  0(AT3). 


We  now  compare  the  throughput  performance  for  OPT 
(i.e.,  exhaustive  search),  Algorithm  1,  and  Algorithm  2.  We 
assume  the  following: 


•  The  path-loss  exponent  is  a  =  3. 

•  The  wireless  channel  is  affected  by  Rayleigh  fading  with 
Rayleigh  parameter  v(S,  D)  =  1. 

•  The  received  power  factor  g{S,  D )  is  given  by  (3). 

•  The  transmit  power  is  Ptx(S)  =  1  for  all  source  nodes  S. 
The  receiver  noise  power  at  destination  D  is  Pnoise(D)  = 
10-4. 


We  now  study  a  stationary  wireless  network  as  shown 
in  Fig.  3,  which  has  a  destination  D  and  K  sources  (Si). 
Assume  that  the  sources  are  located  randomly  in  the  circle 
centered  at  (0,  0)  and  of  radius  r  =  5.  The  destination  is 
located  at  (xn,0).  In  the  following,  we  show  the  through¬ 
put  T  versus  the  SINR  threshold  /3  for  various  network  sizes 
and  topology  configurations.  The  throughput  values  are  av¬ 
eraged  over  100  random  network  instances.  Smaller  values 
of  /3  result  in  higher  throughput  T,  as  expected. 

Consider  a  small  network  with  K  =  10  sources.  First,  let 
Xd  =  0,  i.e.,  the  destination  is  located  at  the  center  of  the 
circular  region  in  which  the  sources  are  distributed.  The  per¬ 
formance  results  are  shown  in  Fig.  4,  which  shows  that  our  2 
heuristic  algorithms  (which  have  polynomial-time  complex¬ 
ity)  perform  almost  identically  to  OPT  (which  is  an  expen¬ 
sive  exhautsive  search).  In  fact,  our  numerical  results  show 
that  Algorithm  1  performs  identically  to  OPT,  and  Algo¬ 
rithm  2  also  performs  identically  to  OPT  with  the  exception 
of  a  few  network  instances. 

Next,  let  Xd  =  10,  i.e.,  the  destination  is  outside  the  circle 
of  radius  r  =  5.  As  shown  in  Fig.  5,  the  throughput  results 
are  lower  than  those  for  the  case  of  Xd  =0.  This  is  because 
the  distances  between  the  sources  and  the  destination  are 
larger  (while  the  receiver  noise  power  at  destination  D  still 
maintains  at  Pnoise(D)  =  10-4),  which  imply  that  the  SINR 
at  the  destination  is  now  reduced.  Again,  our  2  heuristic 
algorithms  perform  almost  identically  to  OPT.  Similar  ob¬ 
servations  (which  are  not  shown  here)  also  hold  for  networks 
of  moderate  size  (e.g.,  I\  =  20  and  K  =  30). 

We  now  consider  a  network  with  K  =  100  nodes.  It 
is  not  feasible  to  apply  OPT,  which  has  high  computational 
complexity,  to  the  network  with  this  large  size.  The  through¬ 
put  results  are  shown  in  Fig.  6  (for  xd  =  0)  and  Fig.  7  (for 
Xd  =  10)  for  our  polynomial-time  Algorithms  1  and  2.  The  2 
algorithms  perform  almost  identically,  with  the  exception  of 
a  few  network  instances  where  Algorithm  1  slightly  outper¬ 
forms  Algorithm  2. 


Fig.  3  A  wireless  network  with  K  sources  (Si)  and 
a  destination  (D). 
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Fig.  4  Throughput  ( T )  vs  SINR  threshold  (/3) 
for  K  =  10  and  xd  =  0. 


Fig.  5  Throughput  (T)  vs  SINR  threshold  (/3) 
for  K  =  10  and  Xd  =  10. 


P 

Fig.  6  Throughput  (T)  vs  SINR  threshold  (/3) 
for  A'  =  100  and  xd  =  0. 


Fig.  7  Throughput  (T)  vs  SINR  threshold  (/3) 
for  I\  =  100  and  xd  =  10. 

Remark  2  One  of  the  algorithms  (OPT,  Algorithm  1,  or 
Algorithm  2)  in  this  section  can  be  used  to  find  a  subset  H 
of  sources  that  transmit  in  a  single  time  slot.  For  the  case 
of  \H\  <  A',  we  can  use  the  algorithm  repeatedly  (with  a 
modified  set  of  available  sources)  to  allow  all  the  sources  to 


transmit  in  multiple  time  slots.  For  example,  suppose  that 
the  total  number  of  sources  is  K  =  6,  and  suppose  that  the 
subset  of  sources  that  transmit  in  each  time  slot  is  found 
to  be  H  =  {aS>2 ,  S5}.  We  can  then  transmit  S2  and  S5  in 
time  slot  1.  We  apply  the  same  algorithm  to  the  reduced 
source  set  {Si,  S3,  S4,  S6}  to  find  the  next  subset  of  sources 
for  transmission  in  time  slot  2,  etc.  In  this  way,  as  a  method 
for  fairness,  all  the  sources  can  transmit  in  multiple  time 
slots.  □ 

7  Summary 

For  the  MAC  problem  considered  in  this  paper,  throughput 
in  a  single  slot  is  maximized  by  the  choice  of  an  appropriate 
transmission  set,  rather  than  via  a  randomized  approach  in 
which  the  sources  transmit  with  some  probability.  This  re¬ 
duces  the  optimization  problem  from  a  continuous  one  to  a 
finite  one.  As  shown  in  our  studies,  the  proposed  heuristic 
algorithms  for  finding  the  transmission  sets  not  only  have  low 
complexity,  but  also  perform  almost  identically  to  the  opti¬ 
mal  sets  obtained  by  expensive  exhaustive  search.  Further, 
our  transmission  method  for  a  single  time  slot  (for  provid¬ 
ing  transmission  opportunities  to  a  subset  of  sources)  can  be 
extended  to  multiple  time  slots  (for  providing  transmission 
opportunities  to  all  sources,  i.e. ,  a  method  for  fairness). 
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